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ABSTRACT 

We present new and stronger evidence for a previously reported relationship between galactic spiral arm pitch 
angle P (a measure of the tightness of spiral structure) and the mass Mbh of a disk galaxy's nuclear supermassive 
black hole (SMBH). We use an improved method to accurately measure the spiral arm pitch angle in disk 
galaxies to generate quantitative data on this morphological feature for 34 galaxies with directly measured black 
hole masses. We find a relation of logiM/Mo) = (8.21 ± 0.16)- (0.062 ± 0.009)P. This method is compared 
with other means of estimating black hole mass to determine its effectiveness and usefulness relative to other 
existing relations. We argue that such a relationship is predicted by leading theories of spiral structure in disk 
galaxies, including the density wave theory. We propose this relationship as a tool for estimating SMBH masses 
in disk galaxies. This tool is potentially superior when compared to other methods for this class of galaxy, and 
has the advantage of being unambiguously measurable from imaging data alone. 
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1. INTRODUCTION 

Since the existence of supermassive black holes (SMBHs) 
as a common, or even ubiquitous, component of galactic 
bulges was first recognized (iKormendv & Richstone 19951 
BarthI 12004 '; Komendy 2004; Magon ian et al.l 11998 ). in- 
creasingly successful attempts have been made to mea- 
sure the mass of these objects. This has enabled as- 
tronomers to discover correlations between the mass of a 
SMBH and its host galaxy's rnass or luminosity ( Kormendv 
1993"; 'Kormendv & Richstone" 1995"; 'MagorrianeLaLl 11998e 
Marconi & Hunt 2003; Hiiring & Rix 2004). A number of 
features of the host galaxy have now been found to correlate 
to the mass of the black hole, giving rise to efforts to study the 
black holes by making measurements of features of the host 
galaxy even where the black hole is undetectable. Although 
many of these correlating features of the host galaxy require 
spectroscopy to measure, on e which does not is the Sersic 
index of the galaxy's bulge dGraham & Driveiil20Q7l). This 
correlation demonstrates the feasibility of estimating SMBH 
masses through imaging data alone. Here we verify and fur- 
ther refine a recently discovered relation between the spiral 
arm pitch angle of a galaxy and the mass of its SMBH, the 
M-P relation (Seigar et al. 2008). 

Our knowledge of SMBH masses in the Universe has grown 
dramatically over the last decade, due primarily to high- 
resolution observations made with the Hubble Space Tele- 
scope (HST). These observations have shown that SMBHs 
reside not only in the cores of active galaxies, as has been be- 
lieved for decades, but also in the centers of quiescent galax- 
ies. Recent works have begun to explore the importance of the 
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nuclear SM BHs in the evolution, or co-evolution, of its host 
galaxy (e.g. Maaonian et al. 1998; Ferrarese & Merritt 2000t 
[Gebhai -dt et al. 2000; Marconi & Hunt 2003; Hiiring & Rixl 
20041; Springel et al. 2005; Hopkins et al. 2007; Rosario et afl 
Ot lCrenshaw et al.ii2010fc iTreuthardt et al.ii2012i) . As a re- 
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suit, any complete theory of galaxy formation h as to produce 
SMBHs in the centers of massive galaxies (e.g. I Silk & ReesI 
1998), and explain the evolution of SMBH mass over time. 

Direct determination of SMBH mass depends on in- 
strumentation which can observe the motion of stars, 
gas, and dust in the immediate vicinity of the black 
hole (Kormendv & Richstone' 1 19951; iMacchetto et alj 119971; 
Macieiewski & Binnev 2001). This process is observation- 
ally expensive, or impossible for distant galaxies. It is nat- 
ural, therefore, that much time and effort has gone into the 
search for an indirect measure which can be used to give 
estimates of central SMBH mass. One remarkable indica- 
tor is the Mbh-c relation, which relates the central SMBH 
mass (Mbh) to the velocity dispersion in the central galactic 
bulge (a) (Ferrarese & Merritt 2000; Gebhardt et al. 2000). 
The Mbh - cr relation has led to considerable success esti- 
mating the SMBH mass for somewhat more distant galax- 
ies, generally those whose cosmological redsh ift (z) is z < 0. 1 
jHeckman et alj2004l;lKauffmann et aU2007l) . Thus the mass 
of the SMBH can be estimated for a wider range of galaxies. 
Nevertheless, as it is necessary to have spectroscopic mea- 
surements to estimate cr for the spheroidal component, it is 
still expensive in terms of telescope time. Additionally, mea- 
suring (7 is more complex for disk galaxies than it is for ellipti- 
cals, since one must account for velocity dispersion associated 
with the motion of disk and bar stars intermingled with the ve- 
locity dispersion of the bulge stars, which is actually the cor- 
relating quantity (Hu 2008). Studies do suggest that the scat- 
ter in the Mbh-c is si gnificantly gre ater for disk galaxies than 
for ellipticals (Giilte kin et al.ll2009h . Although projects such 
as the Sloan Digital Sky Survey have made progress in ac- 
quiring spectra of large numbers of galaxies, there have been 
few ways, up to now, to take advantage of the far larger cat- 
alogs of imaging data available from the public archives for 
the purpose of estimating black hole masses. 

A far greater number of estimates can be made if features 
that are linked to the mass of SMBHs can be measured from 
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imaging data. Several such relations have been explored, in- 
cludi ng those between black hole mass and bulge lum inos- 
ity, (iKormendv & Richstonelll995l:lHaring & Rix"2004, e.g.), 
and S ersic index and nuclear SMBH mass, ([Graham & Driver 
I2007h . among others. More recently, a relation between the 
central SMBH mass and the spiral arm pitch angle of its host 
galaxy has been discovered by examination of 2 7 disk galax- 
ies w ith previously estimated SMBH masses (* Seigar et al.l 
120081) . The pitch angle of the spiral arms of the galaxy (essen- 
tially how tightly the spiral structure in the arms are wound) 
can be measured solely from images of the galaxy. There exist 
sets of images which cover significant field-of-view (FOV) to 
a considerable depth and look-back time for which this tech- 
nique could provide estimates of the SMBH mass in a com- 
plete sample of spiral galaxies at much greater distances than 
has been possible hitherto, up to redshifts z ^ 1 for especially 
deep images and favorable objects, but very likely to a redshift 
of z ~ 0.5 for a significant sample of galaxies. 

We may ask ourselves if such a relationship is expected 
from our understanding of galactic physics. The answer is 
in the affirmative. Firstly, it has been shown empirically that 
there is a link between spiral arm pitch angle and the central 
mass concentration of galaxies (Seigaretal. 2005, 2006). It 
is generally agreed that a "strong correlation between central 
mass concentration an d pitch angle [is] p redicted by modal 
density wave theory," (iGrand et al.l 120 121) the dominant the- 
ory of spiral arm structure. Furthermore such a correlation 
between the size of the central bulge and the tightness of the 
spiral arm winding constitutes essentially the first significant, 
though at that point largely qualitative, observation of extra- 
galactic astronomy, the basis of the Hubble classification. The 
notion that black hole mass depends on the mass of the cen- 
tral galactic bulge is now also widely established, as a result of 
observed correlations of black hole mass with both bulge ve- 
locity dispersion and bulge luminosity. Therefore black hole 
mass and spiral arm pitch angle should each measure the cen- 
tral mass concentration and should therefore correlate with 
each other quite strongly. Although the mechanism by which 
the correlation of black hole mass with central bulge mass is 
maintained is uncertain there are neverth eless highly plausi- 
ble arguments why they should correlate (ISilk & Reeslll998l 
e.E 



IShul (119841) . which deals primarily with the case of den- 
sity waves in the rings of Saturn naturally focuses on what 
can, in the context of galactic astronomy, be termed the bulge- 
dominated case. One has a central body (Saturn) whose mass 
far outweighs that of the material in the disk. The result is rel- 
atively straightforward to deal with theoretically, producing 
density waves in a tractable short wavelength approximation 
(density waves in Saturn's rings typically have pitch angles 
measured in tenths of a degree). Even better the system is 
well understood observationally and there is excellent agree- 
ment between theory and obser vation. We may thus be very 
happy with the result given by IShul (119841) that the tangent 
of the pitch angle (/) of the spiral pattern produced by den- 
sity waves should be proportional to the ratio of the surface 
density of mate rial in the d isk to the total mass of the central 
body. Although Shul (1 1 984 ) focuses mostly on the case of Sat- 
urn the formalism (the theory was, of course, developed in the 
context of galactic spiral arms) is still approximately valid for 
the most bulge-dominated disk galaxies with relatively tightly 
wound spirals (small / and thus with the shortest wavelength 
density waves seen in galactic disks). 

Naturally the case of galactic disks is much more com- 



plicated than that of Saturn's rings, not least because of the 
self-gravity of the disk itself. Nor can we say that there is 
such close agreement between theory and observation in this 
case. Nevertheless there is quite good agreement and some 
success has been achieved in modeling individual galaxies 
with t he density wave theory. An example is Roberts et aLl 
(119751) which shows that the pitch angle in disk galaxies de- 
pends on the ratio of two radii, the half-mass radius, defined 
as the radius within which half the mass of the galaxy's disk 
is contained, and the co-rotation radius, defined as the radius 
at which stars and other material bodies in the disk rotate at 
the same rate as the spiral pattern. Thus, once again, we see 
that the more concentrated the mass of the galaxy is towards 
the center (and thus the smaller is the half-mass radius), the 
tighter will be the spiral pattern. 

It is not hard to show that the result of [Roberts et al.l (119751) 
is compatible with Shu (1984). If one reduces the former's 
Toomre disk model to a very simple bulge (or planet) with a 
thin low mass disk of uniform density and thickness, then the 
half-mass radius shrinks (and thus the pitch angle decreases) 
depending on the ratio of the central mass (the planet or bulge) 
to the surface mass density of the disk, which is the control- 
ling factor in Shu (1984). Thus in bulge dominated galaxies 
the pitch angle depends inversely on the mass of the central 
bulge. Disk dominated galaxies, especially the extreme case 
of bulgeless galaxies, behave similarly, in that their pitch an- 
gle correlates to the relative concentration of mass towards 
the center of the galaxy. Unfortunately relatively little is 
known as yet about the relation between black hole mass and 
galaxy characteristics in disk-dominated galaxies, since rela- 
tively few black hole masses have been directly measured in 
such galaxies. 

Indeed if there is no classical bulge, it is difficult to know 
how to interpret the M-cr or M-bulge luminosity relations at 
all. Clearly, further work will have to be done to understand 
the relation between disk dominated galaxies and their cen- 
tral black holes (if they have one). However there are argu- 
ments which might suggest there should still be a link be- 
tween black hole mass and the mass of the central part of 
the disk. It is widely suspected that bulges are produced by 
mergers, in which the central parts of galaxies become hotter 
and, heated past the point in which they maintain a flattened 
disk profile, adopt a bulge profile. If the central black hole 
mass does indeed correlate to the mass of this merger-created 
bulge, one might speculate that it would have previously (be- 
fore the mergers) correlated to the mass of the central region 
of the disk, out of which material the post-merger bulge was 
presumably formed. But in that case the spiral arm pitch an- 
gle, pre-merger, would have also tended to correlate to the 
mass of the disk's central region (which tends to control the 
value of the half-mass radius). Although it is not possible to 
say anything with certainty at this stage, it may prove that 
the mass-pitch angle relation could work for disk-dominated 
and bulgeless galaxies where other correlations (M-cr, M-L) 
would need to be reinterpreted. In the meantime, we can be 
fairly confident that for galaxies with classical bulges the pitch 
angle of the spiral arms should correlate well to the mass of 
that central bulge. 

The dependence of pitch angle on central mass in the modal 
density wave theory can be understood by analogy with stand- 
ing waves on a string, since the modal density waves them- 
selves constitute a standing- wave pattern. In the case of waves 
on a string one expects the wavelength of the waves oscillating 
between the ends of the string to depend on the speed of prop- 
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agation of the wave. This in turn depends on the ratio p/T 
where p is the density of the string and T is the tension in the 
string, the restoring force producing the wave phenomenon. 
In the modal theory the density of material in the disk cjo plays 
the role of the density of the string and the restoring force or 
tension is produced by the gravitational field of the massive 
central region of the galaxy. It is natural that the wavelength 
of the result spiral pattern should depend on the ratio of these 
two quantities. 

It should be noted that the modal density wave theory is not 
the only theory which attempts to explain spiral arm structure 
in disk galaxies. Rivals in clude the swing amplification model 
of density wave theory (iKormendviri981l: iGerola & SeidenI 
[19781 ISeiden & Gerolalll979l) andthe Manifold theory (see 
below for selected references). It would be fair to say that the 
modal density wave theory is the most widely accepted but 
that each of these has significant support, at least for certain 
types of spiral galaxies. It has even been suggested that dif- 
ferent galaxies (for instance, grand design versus flocculant) 
have different mechanisms explaining their spiral structure. 
Apart from a theory which proposes that spiral arms are the 
result of stochastic star formation operated upon by differ- 
ential rotation all theories agree that there is a link between 
central mass and spiral arm pitch angle. The Toomre density 
wave theory differs from the modal theory primarily in deny- 
ing that spiral arm structure lasts for longer than about a few 
rotational periods in a given galaxy (because in this theory 
there is no longstanding standing wave pattern). Pitch angle 
should vary with time in this theory, but should still obey a 
relation with the size of the central mass. 

The Manifold theory of spiral structure is the most re- 
cently proposed of these theories ( Kaufmann & Contopoulos' 
1996; Harsoula & Kalapotharakos 2009; Athanassoula et al. 
2009blal 120 lOt Athanassoula 2012). This theory describes 



the spiral pattern as being the result of stars formed near the 
ends of a galaxy's bar moving into chaotic, highly eccentric 
orbits which nevertheless cause the stars to move along rela- 
tively narrow tubes known as manifolds. The global pattern 
produced by their motion along these manifolds gives rise to 
the observed spiral arms. The details of this theory are also 
subtle, but it is abundantly clear that the orbits, and therefore 
the manifolds, are controlled by the central mass concentra- 
tion, as with all galactic orbits, and that therefore, once again, 
the pitch angle of the spiral arms should vary with the central 
mass of the galaxy (Athanassoula, private communication). 

Thus, we see that the primary theory for the formation of 
spiral structures, along with its two main competitors, are 
agreed that the mass of a central black hole should correlate 
with the mass of the central core of the galaxy. These theo- 
ries demand that the mass of the central bulge should deter- 
mine the pitch angle of the galaxy's spiral arms. Indeed, it 
is currently difficult to imagine a theory of spiral arm struc- 
ture which did not demand a correlation with the central mass 
concentration, at least indirectly. But the two currently most 
actively pursued theories (modal density waves and manifold) 
both give the central mass a controlling influence on the mech- 
anism which generates the pitch angle of the spiral pattern. It 
is not at all surprising, then, that we should find strong ev- 
idence for such a concentration in actual observations, and 
with a notably low degree of scatter 

Finally, one need look no further than the Hubble sequence 
to see an illustration of the connection between galactic mor- 
phologies and SMBH mass. The SMBH mass-bulge mass 
relation, when combined with the general pattern of larger 



bulges and tighter spiral arms as one moves from Sc to Sa 
in the Hubble sequence, demonstrates, at the very least, an in- 
direct connection between these properties. Our view is that 
spiral arm pitch angle, which appears to be well correlated at 
least with SMBH mass, would be an excellent tool to probe 
the complex of correlated characteristics of spiral galaxies for 
several reasons. First, because it is measurable through imag- 
ing data alone. Second, it can take advantage of the great 
storehouse of publicly accessible archival data available. Fi- 
nally, its measurement is independent of redshift, since log- 
arithmic spirals remain self-similar no matter how they are 
scaled. In short, it may be possible that, for disk galaxies, 
we can gain information on the black hole masses for a sig- 
nificant number of spiral galaxies which previously could not 
have their masses estimated by other means. 

Over the last few decades, it has become widely accepted 
that SMBHs and dark matter play influential, even dominant, 
roles in the evolution of galaxies. As neither black holes 
nor dark matter are directly observable in any part of the 
electromagnetic spectrum, information about them has been 
painstakingly obtained from observations o f their gravita- 
tional effect on baryonic matter. Admittedly, Kormend v et alj 
( 1201 II) suggested that SMBH mass doe s not c orrelate with 
galaxy disks, and iKormendy & Bendol (1201 Ih further sug- 
gests that a galaxy's dark matter halo does not have any direct 
corre lation with the proper ties of the nuclear SMBH. How- 
ever, IVolonteri et d] ( 1201 I I) and 'Booth & Schave' ('2009) both 
provide counter results to th ose of Kormendv et al. (2011) and 
IKormendy & Bended (1201 II) . 

In this paper, we will reexamine the relationship of 
ISeigar et"al] (120081) and expand upon the sample used in that 
study by adding new measurements from other nearby spi- 
ral galaxies and active galactic nuclei (AGN). We will dou- 
ble the number of points used in the previous work, as well 
as update the meth od used in the mea surement of the spiral 
arm pi tch angles in lSeigar et al.l (120081) to that of iDavis et al.l 
( 120121) . We take advantage of SMBH mass data from a vari- 
ety of measurement techniques including direct measurement 
of stellar and gas dynamics in the vicinity of the black hole, 
measurements based on available maser data for several ob- 
jects and reverberation mapping. We also make use of a select 
set of data based on the M-a relation, in an effort to deepen 
our understanding of the extent to which spiral arm structure 
correlates to central mass. 

The structure of this paper is as follows: in Section |2] we 
outline the data we use in this work (including the observa- 
tions and mass determinations) and the techniques used to 
measure morphological features of the observed galaxies. In 
Section [3] we assemble a updated SMBH Mass-Pitch Angle 
Relation using a variety of observational results and compare 
them across different subsamples. We also examine the use 
of Sersic index as a means of estimating SMBH mass from 
galactic morphologies and compare that method with our re- 
sults. In Section II] we discuss the implications of this work 
and the usefulness of this method of SMBH mass estimation. 
Finally, in Section |5] we outline our final assessment of our 
results. 

In this work, where necessary, we assume a cosmology of 
Qa = 0.728, Qh = 0.0455, n,X = 0.1347, and Ho = 70.4 km 
s"' Mpc"'. This corresponds to the maximum likelihood cos- 
mology from the combin ed WMAP-hBAO-hHO results from 
the WMAP 7 data release (iKomatsu et alJl20IlD . 
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2. METHODS 

In this work we investigate the relation between spiral 
arm pitch angle and central black hole mass using galaxies 
selected from a variety of sources with directly measured 
SMBH masses. Candidate galaxy images are selected from 
the available archival data. The sample listed in Table [T] 

includes spiral galaxies that have measurable spiral arm 
pitch angles and measured SMBH masses. For many of these 
galaxies we also measure the Sersic index to compare with the 
results of lGraham & Driveri (12007b . 

2.1. Sample Selection 

The most desirable sample for us to use in this analysis is 
one consisting of spiral galaxies with direct measurements 
of the SMBH mass through examination of either stellar or 
gas dynamics, or both, within the sphere of influence of the 
nuclear SMBH. We have available 10 galaxies with masses 
measured in this way. There are an additional 12 galaxies 
that have upper limits on their measured masses from stel- 
lar or gas dynamics and have no other estimations of their 
masses using other direct techniques as discussed below. We 
do not make use of these limits in constructing our relation, 
but we will later discuss the extent to which these limits are 
compatible with it. Nearly all (7 of 10) of these black hole 
mass measurements are for relatively large black holes, with 
masses Mbh > 1 x 10^ M© in galaxies with relatively tightly 
wound spiral arms, with typically P < 15 degrees. If we notice 
that the mean value for pitch angle in nearby spiral galaxies is 
21.44 degrees (in prep Davis et al. 2013 ). then all but one of 
these measurements, excluding limits, are for galaxies whose 
spiral arms have pitch angle less than this average. This means 
that we are missing the entire right hand side of the distribu- 
tion. That is to say, we have little information on the correla- 
tion for loosely wound spirals. It has been argued that spirals 
with the smallest black holes, i.e. the most loosely wound spi- 
rals, may not fit any co rrelation without considerable scatter 
dKormendv et alj|201 Ih . but substantial improvement in the 
relation would be possible with more data in the region from 
20-30 degrees in pitch angle. There is even more evidence 
that some of these galaxies may not contain black holes at all, 
but merely nuclear star clusters (or black holes within larger 
nuclear star clusters). This, however, is beyond the scope of 
our present discussion. 

Fortunately, other techniques are available which provide a 
considerable number of further data points in order to expand 
our sample. These techniques provide several galaxies with 
smaller mass black holes than are available from the methods 
which can be employed in normal galaxies. We will incor- 
porate data available from maser modeling and reverberation 
mapping. 

The m aser modeling da ta com e s from iLodato & BertinI 
(l2003h: IPastorini et al.1 (l2O07h: llshihara et al.l (120011): 
iKondratko e t alJ (|2006l, 120081); iRodrisuez-Rico et alJ 
(120061): iGre enhill et al. (2003b); Braatz & Gugliucci (2008); 
iGreenhillet al.i (i2003ai) : iKuo et al.i 6201 li) . Observations 
of H2O masers in the vicinity of an active black hole are 
used to obtain circumnuclear disk rotation curves which can 
generate accurate measurements of the mass of SMBHs in 
these galaxies. 

We also use 14 galaxies with mass estimates from rever- 
beration mapping, 12 of which have not been measured us- 
ing other direct methods. For discussions of this method see 
[Peterson e t al,, (2005J; Be ntz et al . (,2009a) . 



Taking these three categories together (direct measurements 
from stellar or gas dynamics, maser modeling and reverbera- 
tion mapping) we have a final sample which includes 10 mea- 
surements using stellar or gas dynamics, 12 using maser mod- 
eling data, and 12 from reverberation mapping. Where the 
three samples overlap, we choose one of them as our preferred 
value (see Table[T]for details). This gives us a final sample of 
34 spiral galaxies with direct measurements of their central 
black hole masses. 

As a check on our work, we will consider a further data set 
of galaxies with more indirect measurements in our discus- 
sion. This includes 4 galaxies with lower limits set by the Ed- 
dington Limit, and 23 (3 also with direct measurem ents) with 
black hole masses estimated by the M-a relation of lFerraresel 
( I2OO2I) . Although there are more recent publications on the 
M-cr relation, we use this source for the purposes of draw- 
ing a comparison with our own previous work and postpone 
a more thorough discussion of the specific relation between 
a and pitch angle to a future work. See Table [T] for full de- 
tails on this extended data set. Also, please note that there is 
overlap between some of the techniques mentioned above. 

In certain cases we have multiple mass estimates or mea- 
surements of galaxies in our data set. Table [Uindicates which 
mass estimates we select for individual galaxies. Here we fa- 
vor direct measurement techniques, such as stellar or gas dy- 
namics, over techniques, such as maser modeling and rever- 
beration mapping, which depend upon the black hole being 
active. 

Whether maser modeling and reverberation mapping 
should be placed in a different category from other direct 
techniques is, of course, highly debatable. Clearly, they be- 
long in a class of techniques which observe signals from ma- 
terial in direct orbit around the black hole itself, rather than 
with techniques such as the M-cr relation, which merely cor- 
relate the mass of the black hole to some feature of the host 
galaxy. But for the purposes of this paper, it is useful to place 
maser modeling and reverberation mapping in a category to- 
gether for two reasons. Firstly, because both methods work 
exclusively for AGN and there has been a recent claim that 
the M-cr relatio n, at least, is d ifferent for AGN than for nor- 
mal galaxies (iParket al.ll2012h . Secondly, because these two 
methods cover a much greater stretch of the sample space than 
the other available direct methods, which tend to have had 
success exclusively for galaxies containing the most massive 
black holes. 

A few exceptions to favoring direct stellar/gas based mea- 
surements exist. In cases where these are available we have 
chosen the most recent and reliable mass estimates available 
for the galaxy from among a variety of measurement tech- 
niques. Where these only produce an upper limit but not an 
estimated mass, we have chosen, where available, mass de- 
terminations from another method which provides a measure- 
ment of the mass with errors instead of just a limit. These dif- 
ferences are noted in Table [T] In Table [T] we include multiple 
entries for each object which has multiple different measure- 
ment types available in the literature. Those that we prefer are 
labeled in the final column of the table. 

Besides these general choices other exceptions ar e also 



made. NGC 5055 has a direct mass estimate in GUI te kin et al 
(J2009) and iBlais-Ouellette et al.1 (12004 . iGultekin et aT 



( I2009.) suggests that the modeling used in the mass deter- 
mination is very uncertain. In this case we have chosen to 
fall back on the mass of the black hole derived from the M-cr 
relation in .Ferrarese , (2(302). In the case of NGC 4395, we 
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TABLE 1 
Sample Data 



Galaxy 


P (deg.) 


Image Source 


Filter 


log(MBH/A^0) 


Measurement Type 


Source 


PrefeiTed 


3c 120 


10.7±1.4 


Danish 1.54 (NED") 


R 


7.72ti 


Reverberation Mapping 


1,2 


Y 


Ai-k 120 


5.4 ±0.6 


HSTACS(Bentz*) 


F550M 


8.15tll 


Reverberation Mapping 


1 


Y 


Circinus 


26.7 ±5.0 


HST-WFPC2 (NED°) 


F814W 


£ 24+0.07 
° -0 08 


Maser ModeHng 


3 


Y 


IC342 


23.2 ±2.8 


VLA (NED") 


21cm 


< 5.70 


Stars/Gas 


4 












6 32+0"** 


M-cr 


5 


Y 


IC 2560 


16.3 ±6.4 


2.5 m du Pont (CGS'') 


B 


f. (;4+0.30 
7.59|i 


Maser ModeHng 


6 


Y 


M31 


8.5 ±1.3 


GALEX (NED") 


NUV 


Stars/Gas 


7 


Y 










M-cr 


5 




M33 


34.5 ±8.6 


Spitzer IRAC (NED") 


IRAC 3.6 ^m 


<3.48 


Stars/Gas 


8 












4 94+0.08 

7 nnSl'M 
'•'"-0.16 


M-cr 


5 


Y 


Mrk 590 


8.5 ±3.3 


HSTACS(Bentz*) 


F550M 


Reverberation Mapping 


1 


Y 


Mrk79 


13.2±3.6 


Lick Im (NED") 


V 


Reverberation Mapping 


1 


Y 


Mrk 817 


9.9 ±4.2 


HSTACS(Bentz*) 


F550M 


131 

<: 04+8:08 


Reverberation Mapping 


2 


Y 


Milky Way 


22.5 ±2.5 


Leiden/ Argentine/Bonn (LAB) Survey (NED") 


21cm 


Stars/Gas 


9 


Y 










M-cr 


5 




NGC 0253 


17.9 ±2.0 


2MASS 1.3m (NED") 


Ks 


yQj+0.30 
7 99+}):08 

7.85_(,29 


Maser ModeHng 


10 


Y 


NGC 0753 


13.2±0.6 


INT 2.5m (NED") 


B 


M-cr 


5 


Y 


NGC 1068 


20.6 ±4.5 


UKSchmidt (NED") 


468 nm 


Maser Modeling 


11 


Y 


NGC 1300 


10.3 ±1.8 


2.5 m du Pont (CGS'') 


B 


Stars/Gas 


12 


Y 


NGC 1353 


13.7±2.3 


2MASS 1.3m (NED") 


K, 


6.68^:«| 


M-o- 


5 


Y 


NGC 1357 


11.8±4.8 


2.5 m du Pont (CGSO 


B 


7 99+0.08 
' "^^-0.09 


M-o- 


5 


Y 


NGC 1417 


12.9±4.1 


2.5 m du Pont (CGS") 


V 


7.62+''-''* 
g g8+<):02 

i7+}):iT 
»-l'-0.15 


M-o- 


5 


Y 


NGC 2273 


17.5 ±7.2 


KPNO 2.1m (NED") 


K 


Maser ModeHng 


13 


Y 


NGC 2639 


12.9±1.2 


HST (NED") 


F606W 


M-o- 


5 


Y 


NGC 2742 


32.5 ±7.9 


CAHA 2.2m (NED") 


1.25 Atm 


J g2+0''8 

g 00 
g9g+}j:ig 

'■"°-0.02 


M-o- 


5 


Y 


NGC 2841 


7.1 ±1.5 


Spitzer IRAC (NED") 


IRAC 3.6 Aim 


M-o- 


5 


Y 


NGC 2903 


15.1 ±3.0 


Pal 60inch (NED") 


440 nm 


M-o 


5 


Y 


NGC 2960 


7.5 ±1.7 


Palomar 48-inch Schmidt (NED") 


645 nm 


Maser Modeling 


13 


Y 


NGC 2998 


14.5 ±9.4 


KPNO 2.1m (NED") 


656.3 nm 


7 ns+0.08 


M-cr 


5 


Y 


NGC 3031 


15.4±8.6 


Spitzer IRAC (NED") 


IRAC 5.8 Atm 


7.91+"" 
1 g5+<):o9 

6 \a^'-^^ 

°-'"-0.09 


Stars/Gas 


14 


Y 


NGC 3145 


7.2±1.3 


2.5 m du Pont (CGS") 


B 


M-o- 


5 


Y 


NGC 3198 


30.0 ±6.7 


KPNO 2.1m CFIM (NED") 


700 nm 


M-o 


5 


Y 


NGC 3223 


10.9 ±2.2 


2.5 m du Pont (CGS") 


B 


7 01 +0.09 


M-o- 


5 


Y 


NGC 3227 


12.9 ±9.0 


JKT (NED") 


Ha 


7 ,0+0. 18 
'■""-0.24 


Stars/Gas 


15 


Y 










Reverberation Mapping 


16 




NGC 3310 


22.7 ±9.1 


HST WFPC2 (NED") 


F814W 


<7.62 


Stars/Gas 


17 


Y 


NGC 3351 


11.1±1.8 


CTIO 4.0m (NED") 


B 


<6.96 


Stars/Gas 


18 


Y 


NGC 3367 


36.8 ±5.3 


CAN Martir 2.12m (NED") 


I 


>5.20 


Eddington 


19 


Y 


NGC 3368 


14.0 ±1.4 


VATT Lennon 1.8m (NED") 


R 


6.90!«:?g 


Stars/Gas 


17,20 


Y 



have decided to use the more recent reverberation mapping 
data and corresponding mass estimate over the mass estimate 
based on combining up per and lower limits set on the mass in 
iFihppenko &H0I (I2003I) . 

2.2. Imaging Data 

The majority of images used for measuring pitch angles 
came from the NASA/IPAC Extragalactic Database (NED)B 
This resulted in the implementation of a wide range of wave- 
length images; anywhere from far-ultraviolet (FUV) to 21 cm 
H emission. See Table[T]for details on each individual galaxy. 
Despite this broad range of wavelength imaging, recent results 
show that galactic pitch angle measurements are i ndependent 
of the wavelength of the image (ISeigar et al.l2 0Q6. Davis et al. 
[2012), or at least not strongly dependent (JGrosbol & Patsis 
[1998) , especially in the UV to near-IR wavelength regimes. 
Although this correlation has not been tested in the mid-IR to 
radio regimes, we can assume that it still applies for several 
reasons. For the mid- IR regime, i.e., where we have used 
Spitzer/IRAC 5.8-/im and 8.0-/im imaging, star formation is 

'^ Ihttp: //ned. ipac. caltech.edu/| 



being traced, and so this should give a similar pitch angle to a 
B-band image. In two cases we resort to 21 -cm data to mea- 
sure pitch angles. While there is currently no empirical ev- 
idence correlating 21-cm pitch angles to optical or near-IR 
pitch angles, there is no reason to suggest that the 21-cm data 
is not influenced in a similar way by the underlying density 
wave. 

When given the option, the imaging with the best resolution 
was used (typically B-band images). According to iThomJevI 
( 1996), a spiral that appears flocculant in the B-band may ap- 
pear to have a weak grand design spiral in the near-infrared 
(NIR). In these case s, near-IR imaging was investigated (typ- 
ically from 2MASS. Ijarrett et all (l200C)h . or Spitzer images). 

Some galaxies in our sample had pitch angles measured 
in Davis et al. (2012), where the method of measuring pitch 
angle we use here is described. S ome of the s e prev iously 
measured pitch angles, reported in iDavis et aTl (12012b . used 
high-qualit y imaging frorn the Carnegie-Irvine Galaxy Sur- 
vey (CGSO; iHo et al.ll2011i) . providing a desirable set of input 
imaging for our two-dimensional Fast Fourier transform sof t- 

^ |http: //cgs ■ obs . carnegiescience .edu/l 









Benier et al. 
















TABLE 1 
















Sample 










Galaxy 


P (deg.) 


Image Source 


Filter 


log(MBH/MQ) 


Measurement Type 


Source 


PrefeiTed 


NGC 3393 


13.1 ±2.5 


CTIO 0.9m (NED") 


B 


7.52:!]« 


Maser Modeling 


21,22 


Y 


NGC3516 


10.6±4.3 


HST WFPC2 (NED") 


500.7 nm 


7.6i:»i« 


Reverberation Mapping 


1 


Y 


NGC 3621 


12.7±1.2 


2.5 m du Pont (CGS') 


B 


>3.60 


Eddington 


23 


Y 


NGC 3783 


10.5±4.8 


LCO 2.5m (NED") 


K 


7.45tli 


Reverberation Mapping 


1 


Y 


NGC 3938 


22.4 ±7.2 


KPNO 2.1m CFIM (NED") 


B 


>4.26 


Eddington 


19 


Y 


NGC 3982 


14.0±0.4 


MaunaKea2.24m (NED") 


R 


<7.93 


Stars/Gas 


18 


Y 


NGC 3992 


6.2 ±6.1 


MaunaKea2.24m (NED") 


B 


<7.78 


Stars/Gas 


18 


Y 


NGC 4041 


23.3 ±8.2 


Palomar 48-inch Schmidt (NED") 


645 nm 


<7.33 


Stars/Gas 


24 


Y 


NGC 4051 


29.1 ±4.9 


MaunaKea2.24m (NED") 


B 


^24+0.12 
"•^^ 0.16 


Reverberation Mapping 


32 


Y 


NGC 4062 


12.4±1.4 


1.8m Perkins (NED") 


B 


"•"•'-0.09 


M-cr 


5 


Y 


NGC 4151 


11.8±1.8 


VLA (NED") 


21 cm 


7 6fi+<'"5 
/.DO 05 


Stars/Gas 


25 


Y 








21 cm 


7.64+"" 
7 cn+o:m 


Reverberation Mapping 


16 




NGC 4258 


7.7±4.2 


Spitzer IRAC (NED") 


IRAC 8.0 iim 


Stars/Gas 


17 


Y 










Maser Modehng 


26,27 












7.48+008 


M-cr 


5 




NGC 4303 


13.5 ±4.6 


1.3m McGraw-Hill (NED") 


B 


6.92!»:» 


Stars/Gas 


17 


Y 


NGC 4321 


21.8±3.6 


KP 2.1m CFIM (NED") 


R 


< 7.46 


Stars/Gas 


18 










R 


g47+o.08 
6 9S+^-^? 
556^-^i 

•'•-'"-O 16 


M-cr 


5 


Y 


NGC 4388 


26.2 ±8.2 


KPNO 2.3m (NED") 


K, 


Maser Modeling 


13 


Y 


NGC 4395 


35.2 ±6.8 


GALEX (NED") 


FUV 


Reverberation Mapping 


28 


Y 


NGC 4450 


9.1±3.1 


KP 2.1m CFIM (NED") 


B 


<8.07 


Stars/Gas 


18 


Y 


NGC 4501 


12.7±1.4 


KP9 t2ka (NED") 


R 


<7.98 


Stars/Gas 


18 


Y 


NGC 4536 


14.8±7.9 


KP 2.1m CFIM (NED") 


B 


>3.68 


Eddington 


19 


Y 


NGC 4548 


25.6 ±6.6 


JKT (NED") 


B 


<7.55 


Stars/Gas 


18 


Y 


NGC 4593 


20.2 ±2.7 


2.5 m du Pont (CGS'') 


I 


6.97t» 


Reverberation Mapping 


33 


Y 


NGC 4800 


21.5±3.2 


KP9 t2ka (NED") 


R 


< 7:53 


Stars/Gas 


18 


Y 


NGC 5033 


16.5 ±5.6 


KP 2.1 CFIM (NED") 


B 


7 24.+0.08 


M-cr 


5 


Y 


NGC 5055 


14.9 ±6.9 


Spitzer IRAC (NED") 


IRAC 5.8 fim 


ggO+«.08 
7 03+^-'^ 
7:80|j|| 


M-cr 


5 


Y 


NGC 5495 


27.8±1.2 


UK 48-inch Schmidt (NED") 


468 nm 


Maser Modehng 


21 


Y 


NGC 5548 


15.0±2.5 


HST (NED") 


F606W 


Reverberation Mapping 


34,35 


Y 


NGC 6323 


11.8±3.4 


Palomar 48-inch Schmidt (NED") 


645 nm 


697+0.00 
6 77+^'^ 

7:5o|i 


Maser Modeling 


13 


Y 


NGC 6926 


17.5±5.5 


2MASS (NED") 


Ks 


Maser Modehng 


29 


Y 


NGC 7331 


22.2 ±4.2 


Spitzer IRAC (NED") 


IRAC 3.6 fim 


M-cr 


5 


Y 


NGC 7469 


28.5 ±4.3 


HST NIC2 (NED") 


FllOW 


<7.73 


Stars/Gas 


16 












7.06+*" 
'•'^-0.18 


Reverberation Mapping 


1 


Y 


NGC 7582 


14.7 ±7.4 


ESQ Im Schmidt (NED") 


R 


Stars/Gas 


30 


Y 


NGC 7606 


11.3±1.2 


2.5 m du Pont (CGS') 


V 


7 27+0-08 

(. Qf:+S:§8 
"•^°-0.26 


M-cr 


5 


Y 


UGC 3789 


10.5±4.8 


Palomar 48-inch Schmidt (NED") 


645 nm 


Maser Modeling 


31 


Y 



Note. — Column 1 Galaxy name. Column 2 Pitch Angle (degrees), The Pitch Angle for the Milky Way comes fromlLevine et al. '200<5?). Column 3 Image Source 
a Images taken from NASA Extragalactic Database (NED), b Images provided by Misty Bentz. c The pitch angles for Cai^negie-Irvine Galaxy Survey (CGS) galaxies 
come from Davis et al. i2012). Column 4 Image filter. Column 5 log(Mass) (Mbh/A^o). Column 6 Measurement type. Column 7 mass measurement source. Column 
8 PrefeiTed measurement. References: (1) Peterson et al. (2004), (2) Bentz et al. 1 2009b), (3) Greenhill et al. (2003a), (4) Bokeretal. (1999), (5) FeiTarese (2002), 
(6) Lshihara et al. (2001), (7) Bender etal. 12005), (8) Meiritt et al. (2001), (9) Gillessen et al. (2009), (10) Rodriguez-Rico et al. (2006), (11) Lodato & BeiTin (2003), 
(12) Atkinson et al. ( 2005), (13) Kuo et al. (201 1), (14) Devereux et al. ( 2003), (15) Davies et al. ( 2006), (16) Hicks & Malkan ( 2008), (17) Pastorini et al. ( 2007), (18) 
gajzietal. (2002), (19) Satvapal et al. (2008), (20) Nowak et al. (2010b), (21) Kondratko et al. i2006), (22) Kondratko et al. (2008), (23) Satvapal et al. (2007), (24) 

(29)lGreenhill et all l2003bl) . (30) 



Marconi et all |2003|) , (25) Onken et al. ( 2007), (26) Herrnstein et al. ( 2005), (27) Mivos hi et alj |1995|\ ( 281lPeterson et al.l |2005|). (29)lGreenhill et al 
Wold etalj 120061) . (31) Braatz & Gugliucci ( 2008) (32) Denney et al. (2010) (33) Dennev eTaTj 120061) (34) IBentz et alj 120071) (351 IBentz et alj 12009^ 



ware, named 2DFFT. Images highlighting the sharpest detail 
in spiral arm struc ture were primaril y selected (typically B- 
band images). See lDavis et alJ ( 120121) for details on these im- 
ages. 

Finally, for some galaxies with active galactic nuclei in our 
sample, we have consulted high-resolution HST ACS F550M 
images of r e verbera tion-mapped AGN host galaxies, used in 
iBentz et al] (l2009ah . Many of these images did not reveal 
any spiral structure because of the bright nucleus or the small 
FOV. In these cases we resorted to ground-based, wider FOV 
images to measure pitch angles. In the end, th ree galaxies 
had pitc h angles measured using HST data from lBentz etalJ 
(l2009ah . 



2.3. Measuring Pitch Angles 

We use the method described in i Davis et al] (120121) to ac- 
curately measure the pitch angles of the 34 galaxies with di- 
rect mass measurements that comprise our sample, plus an 
additional 33 galaxies with mass limits or M-cr estimates used 
in our extended data set. This technique is an extensio n 
of the method described in Sarai va Schroeder et al.l (11994b . 
which utilizes a 2DFFT algorithm to measure the pitch an- 
gle between a user defined inner and outer radius on a de- 
proiected image of a sp iral galaxy. For more details see 
iPuerari & Dottoril (119921) : iPuerari et^D (120001) . Galaxies are 
deprojected by assuming that the disk galaxy, when face-on, 
will have circular isophotes. Although Ryden (2004) has 
shown that disk gala xies do have an in trinsic ellipticity, it is 
relatively small, and iDavis et alJ (120121) has shown that small 
errors in the measured axial ratio of galaxies do not affect the 



SMBH Mass - Pitch Angle Relation 



measured spiral arm pitch angle. 

The extension to this method in iDavis et aTl (120 12h elimi- 
nates the user defined inner radius in favor of measuring the 
pitch angle over all possible inner radii, thus allowing the user 
to examine the results of the Fourier analysis for long regions 
over which changing the inner radius of the transformed re- 
gion of the data does not affect the measured pitch angle. This 
provides a far more accurate measurement of the pitch angle 
of the galaxy than single measurement techniques utilizing in- 
dividual inner radii, as well as providing us a means of exam- 
ining the consistency of the logarithmic structure of the spiral 
arms. For further details on this technique see iDavis et aP 

There are many advantages to using this method. First, it 
helps us eliminate a great deal of uncertainty involved in mea- 
suring the pitch angles of galaxies. Instead of simply assum- 
ing that the galaxy's spiral is logarithmic, it provides some 
check on the extent to which that is true by varying the re- 
gion over which the Fourier analysis is performed. If a depar- 
ture from logarithmic behavior is found (most commonly in 
a change of the spiral arm structure in the very outer regions 
of the galaxy) then the user can select a region over which the 
behavior is consistent and use the pitch angle associated with 
that region. In our case this means we are able to focus on 
the inner region where we expect that the behavior of the spi- 
ral arm structure will be more strongly affected by the mass 
of the central concentration. Additionally, both iSeipar et alj 
(12006) and Davis et al. (2012) established the consistency be- 
tween B-band and NIR-band pitch angles, demonstrating that 
pitch angle does not depend measurably on the band chosen 
for imaging. The results of IDavis et al.l (120 12h demonstrate 
that the pitch angles of the galaxies in that sample appear to be 
generally independent of the band of the observations, at least 
within the uncertainties of the reported pitch angles. Only one 
group have so far reported such an effect (Grosbol & Patsis| 
Il998h and they agree that the amount of the variation (which 
is visible in only three of their studied galaxies) is no more 
than seven degrees between K band and B band images. Thus, 
with this method one may measure the spiral arm pitch angle 
for nearly all late-type galaxies, limited by little more than the 
requirement, which should be random and unbiased, that the 
galaxy is not close to edge-on to our line-of-sight. 

One obvious variation in the appearance of spiral arm 
galaxies is the distinction between flocculant and grand design 
spirals. Flocculant spirals lack the regularity in arm segments 
seen in other spirals. Because of this lack of long stretches of 
continuous arms, measuring a pitch angle for them is not so 
straightforward, at least using manual methods, as it is in the 
case of grand design spirals. However, our method permits 
the user to establish whether there is, nevertheless, a consis- 
tent spiral pattern, with measurable pitch angle and number 
of arms, over the whole disk. Only in a handful of cases in 
this sample does a flocculant galaxy present particular dif- 
ficulties to our method, as evidenced by a larger than usual 
error quoted. 

It has been argued (for instance bv lSeiden & Geroial(ll982l) ) 
that flocculant spiral arms may be produced by a different 
physical process than grand design spiral arms (regions of 
self -propagating star-formation acted upon by differential ro- 
tation in the former case, spiral density waves in the latter). 
It is also noteworthy that galaxies which are flocculant in 
the B-band may show a grand design pattern in the R-band 
(lThomlevll996h due to the potentially different origins of the 
spiral structure in flocculant galaxies and the old stellar pop- 



ulation in the redder b ands tracing the spiral density waves 
dSeigar & James|[T998l) . It appears that though there may be 
a large difference in detailed appearance (thus in one band, 
the image seems flocculant, but not in another) there are only 
small differences in measurable structure, the spiral arm pitch 
angle. This is certainly suggestive of a similar underlying 
cause for the different varieties of spiral structure. In our case, 
where we did find large errors in the measurement of floccu- 
lant spirals we have preferred to use near-infrared images for 
such galaxies, which seems to reduce measurement errors no- 
ticeably. 

2.4. Sersic Index 



Sersic index (ISersid Il963h has been proposed as an- 
other o bservable feature which correlates with SMBH mass 
d Graham & Driver 2007). Like Pitch Angle, it can be mea- 
sured using only imaging data. It may be that a combina- 
tion of these two approaches. Pitch Angles for spiral galaxies 
and Sersic index for ellipticals, lenticulars and edge-on spirals 
will enable observers to estimate the SMBH mass function 
in normal galaxies out to considerable distances. Although 
Sersic indices can be measured for face-on spirals, doing so 
involves a complex process of disentangling bulge from disk 
and bar components of the galaxy. It is likely that measuring 
the pitch angle of such galaxies will be easier and more ac- 
curate. However, any such campaign will certainly demand 
some analysis of how well the two approaches agree in their 
estimate of central black hole mass. In order to make a com- 
pa rison between our pitch a ngle-derived masses and the work 
of iGraham & Driverl (12007 ). we must calculate the Sersic in- 
dex for the galaxies in our sample. We have done so for 
3 1 galaxies, and a further 4 have been taken from the liter- 
ature! Graham & Driver 2007; Kent et al. 1991; Nowak et all 
l2010ah . We have excluded AGN from the sample of galaxies 
with measured Sersic index due to the difficulties presented 
by the bright nucleus. 

The Sersic profile relates how the brightness of a galaxy 
falls off with distance from the center It is of the form; 



/(/?) = /,exp 



-i„[(«/fi.)"-i] 



(1) 



where R is the radius of the isophote, R^ is the radius which 
encloses half of the light of the galaxy, le is the intensity at 
this radius, b„ is fitted with the function: 

^=1.9992«- 0.3271 (2) 

( iGraham & Driverl 12005b . and lastly, n is the Sersic index. 
The Sersic index is also a measure of the concentration of 
the galaxy, defined as the amount of light enclosed by some 
fraction (usually taken to be around a third) of the effective ra- 
dius divided by the amount of light enclosed by the effective 
radius, which by definition is half the light of the bulge. 

We measure the Sersic index by fitting the surface bright- 
ness contours produced by IRA F, Image Reduction and Anal - 
ysis Facility, routine Ellipse (Tody 1986. lJedrzejew ski 1987). 
A bulge/disk decomposition is performed on this data and a 
Sersic profile (plus exponential disk for spirals) is fit to the 
resultant data. 

It was shown in IGraham et al.l (120011) that light concen- 
tration correlates strongly with black hole mass. Later, 
Graham & Driver (2007) found a log relation between Sersic 
index and black hole mass. This correlation is of the form 



/og(Mbh) = (7.98 ± 0.09) + (3 .70 ± 0.46)log(n/3) 
-(3.10±0.84)[;og(n/3)]^. 



(3) 



Benier et al. 



TABLE 2 
Comparison to Sersic Index 



10" - log(Mj„/MQ) = (8.21±0.16) - (0.068±0.009)P - 



Galaxy Moq^hology Sersic Index logiMgn/MQ ) Source 



IC 2560 
M31 


(R')SB(r)b 
SA(s)b 


i.i9:p* 


6.994^-0« 
8.08ltij 




Milky Way 


S? 


1 12+0.22 


7 055+0.081 


2 


NGC 0253 


SAB(s)c 


, ,7+0.23 


f, 005+0.039 




NGC 0753 


SAB(rs)bc 


l-5»-o:26 


7 ,00+0.093 
'■1^^-0.119 




NGC 1357 


SA(s)ab 


1 /Q+fl.30 
^■™-0.25 


7.137+0.083 

7.17g+S:(')§I 

8 128 

7 577+0:199 

'■3' '-0.203 
im+0.40l 

7 ,qo+0093 
'■^^'J-0.119 




NGC 1417 


SAB(rs)b 


1.56+0.31 

3.28 „„ 
9 97+0.45 
^■^'-0.38 




NGC 2903 


SAB(rs)bc 




NGC 2998 


SAB(rs)c 




NGC 3031 


SA(s)ab 


3.23-^it 


3 


NGC 3145 


SB(rs)bc 


l-58!g:i 




NGC 3198 


SB(rs)c 


1.33!«:| 


7 060+<'060 
'■"°"-0.091 




NGC 3227 


SAB(s)a pec 


2 52+0-50 
1 sQ+o:38 

'■'^^-0.32 


7.71g+0.235 

7.362+o:i44 

6.914+0:023 
7 ^oo+o:™ 




NGC 3310 


SAB(r)bc pec 




NGC 3351 
NGC 3367 
NGC 3368 


SB(r)b 

SB(rs)c 

SAB(rs)ab 


2 40+0-**^ 
o:98+j| 

^■-^ -0.00 


4 


NGC 3621 


SA(s)d 


1.89+0 3^ 
J 45+0:25 

2J4+8:M 

1 40+^^-^^ 


7 362+<''l'*'' 




NGC 3938 


SA(s)c 


7.121+0.076 

7 5o4+o:i8D 
7^095^:*! 




NGC 3982 
NGC 3992 


SAB(r)b: 
SB(rs)bc 




NGC 4041 


SA(rs)bc 


74+015 


6-876^:S^i 




NGC 4051 


SAB(rs)bc 


2 ig+0.44 
2.04+0:34 


7.527+0.182 
7 447+S:i{)^ 

6-877^:1 




NGC 4258 


SAB(s)bc 




NGC 4303 


SAB(rs)bc 




NGC 4321 
NGC 4450 


SAB(s)bc 
SA(s)ab 


1.34+fl.27 

2.28+0'*^ 
T 07+0:17 

'■="^-0.26 


7.345+0.136 

7.065+^:^^''' 
7 533+0:199 
7.577+8:8§? 
7 190+0:053 

'■^^"-0.119 




NGC 4501 
NGC 4536 


SA(rs)b 
SAB(rs)bc 




NGC 4548 


SB(rs)b 




NGC 4593 


(R)SB(rs)b 




7.651!°;^] 




NGC 5033 
NGC 5495 
NGC 6926 
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Note. — Column 1 Galaxy name. Column 2 Galaxy Moiphology taken from 
NED. Column 3 Sersic index. Column 4 log(Mass ) (Mbh/A/q ) derived from each 
galaxy's measured Sersic index using the relation of lGraham & Driver 1 2007 ). Col - 
u mn 5 Source of Sersic inde x m easurement. (1) is this work, (2). Kent et al, 1 199 II) . 
OI IGraham & Driveil<2007l) . (4) INowaket"ai])2010at . 

Thus the Sersic index provides an estimate of SMBH masses 
through images of galactic bulges. We will make comparisons 
betwe en our technique and the results of iGraham & Driver! 
(I2007h in Section 13.71 The Sersic indices for several galax- 
ies in our sample, along with corresponding mass estimates 
from the relation of lGraham & Driven (120071) are included in 
Table |2] 

3. RESULTS 

3.1. An Updated SMBH Mass - Pitch Angle Relation 

As discussed above, we chose our sample to include spiral 
galaxies whose central black hole masses have been measured 
using a direct technique. We define a direct technique to be 
one which measures the motions and positions of material in 
orbit around the black hole or directly within its sphere of 
influence. This definition encompasses quite a few different 
methods. Not all are equally accurate or reliable, but for our 
purposes they all have the important distinction that they do 
extract information from signals emitted by material in the 
direct gravitational influence of the black hole. 
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Fig. 1. — Black hole mass versus pitch angle for all spiral galaxies with 
directly measured black hole masses available. The best linear fit to this data 
is illustrated and gives the following relation /og(MBH/'WQ) = (8.21itO. 16)— 
(0.062 ± 0.009)P. The fit has a reduced x^ = 4.68 with a scatter of 0.38 
dex. Black hole masses measured using stellar and gas dynamics techniques 
are labeled with black X points (10 points), reverberation mapping masses 
with red triangles (12 points), and maser measurements with blue squares 
(12 points). 

Combining these three samples, stellar and gas dynamics, 
maser modeling, and reverberation mapping, we have a total 
sample of 34 objects. Fitting these data points (Figure [T]!, we 
find an updated SMBH M-P relation of 

/o§(Mbh/M0) = (8.21±O.16)-(O.O62±O.OO9)P. (4) 

where P is the absolute value of the measured pitch angle of 
the galaxy in degrees. Please note that the sign of a pitch an- 
gle measurement merely represents the chirality of the galaxy 
based upon the user's line-of-sight and is unimportant. This 
fit has a x^ = 4.68 and a scatter of 0.38 dex. This res ult is 
consistent with the previous result of lSeigar et al.l (120081) . It is 
encouraging that, with significantly more data available, the 
scatter has remained unchanged. This value is less than the 
scatter o f spiral galaxie s abou t the M-ct relation of ^ 0.56 
dex from lGiiltekin et akl (120091) . A Pearson Rank Correlation 
Coefficient test produces a coefficient of -0.81, a strong anti- 
correlation. This result has a significance of 99.7% a 3-cr 
result. 

As discussed above, the scatter in this relation is com- 
parable to other relationships. A possible reason for the 
reduced scatte r in this relationship, when compared with 
iGultekin et al.l ([2009), is that the measurement of a in spiral 
galaxies requires that one distinguish the contributions of the 
galactic bulge from other stellar components such as the disk 
or bar Where the galaxy has an active nucleus, the region 
over which a is measured is usually obscured and a proxy 
(some spectral line or lines in the AGN spectrum) must be 
used. Many of the galaxies in this sample have AGNs in their 
nuclei. 

In the case of bulge luminosity the presence of an AGN 
can sometimes be overcome, but when the galaxy is a disk 
galaxy one must undertake decomposition of the light from 
the galaxy, fitting multiple components to the luminosity pro- 
file (bulge, disk and bar are the main components). Once 
again, the presence of unresolved nuclear flux (e.g., AGN or 
LINER or even nuclear star cluster) will create difficulties in 
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determining bulge luminosity or Sersic index. Some algo- 
rithms have been used to measure large numbers of galactic 
bulge luminosities for the purpos es of bulge mas s estimates 
or measurements of Sersic index. (Ide Souza et al.l l2004. e.g.). 
Automated codes may have difficulties disentangling the un- 
resolved nucleus and the bulge, especially with ground-based 
se eing of ~ 1" for an y galaxies that aren't nearby. The method 
of lDavis et al.l (120121) is unaffected by this issue. In the case of 
measurements of spiral arm pitch angle, we examine a com- 
ponent which is unambiguous in spiral galaxies in much the 
same way t hat a is unambigu ous in elliptical galaxies. The 
method of (iDavis et al J 120121) requires minimal image pro- 
cessing to yield a measurement of the pitch angle, requiring 
only a deprojection and cropping of the image and the loca- 
tions of the central pixel and galactic edge. 

For the purposes of estimating black hole masses on the 
basis of pitch angle measurements of galactic spiral arms, we 
propose to use this fit based on direct measurements of SMBH 
mass. 

3.2. Comparison with the previous M-P Relation 

In our previous work ("Seigar et al.l2008l ). black hole masses 
derived from measurements of a, th e velocity dispersion in 
the core bulge region, all reported in iFerraresel (£002 ), were 
used. This was largely because of the relative scarcity of di- 
rect measurements of black hole masses in spiral galaxies. 
Since then the available number of direct measurements has 
increased by over a factor of two. We therefore exclude any 
cr-derived masses from the correlation fit developed here, as 
it is not a direct measurement of black hole mass. Indeed, 
the correlation between a and black hole mass which we use 
is based upon a fit using some of the objects we include in 
our sample. Later in this section we will include those same 
cr-derived values in a fit which we derive primarily to check 
how closely our resul ts agree with the results reported in that 
earlier paper lSeigaret al. (2008). Additionally, when we com- 
pare our results with these indirect values it provides a useful 
check on the overall validity of our correlation. This is espe- 
cially true because the number of black hole masses measured 
by direct techniques which are near or below a million solar 
masses is very small. Including some a values does provide 
a useful check on the slope of the fit by providing some extra 
evidence at the low mass end of the graph. In any case the 
value of our fit, including the cr-derived values, agrees very 
closely with the one reported in the earlier paper 

It is important to note that while the updated fit reported 
here differs from the one in Seigar et al. (2008), though only 
to a minor degree, the discrepancy is not caused by the re- 
measurement of the pitch angles by our improved technique, 
or by revisions of the black hole mass measurements used pre- 
viously. Both of these changes were quite minor in any case. 
More importantly, it is not a new trend indicated by the new 
black hole mass measurements which were unavailable be- 
fore. Rather, the difference is purely because we are in a posi- 
tion to dispose of the use of masses derived by indirect meth- 
ods. It should be understood, and it will be shown below, that 
the results of this paper agree remarkably closely with those 
previous results when we include those indirect masses used 
previously. It is only by removing these that we differ at all 
from the earlier result. 

As we will argue in Section ID we see the relation between 
pitch angle and black hole mass as a natural result of density 
wave theory, which demands that the wavelength of the spi- 
ral density waves should depend directly on the mass of the 
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Fig. 2. — Black hole mass versus pitch angle for all spiral galaxies with 
directly measured black hole masses available. Black hole masses measured 
using stellar and gas dynamics techniques (10 points) are labeled with black 
X points, reverberation mapping masses with red triangles (14 points), and 
maser measurements with blue squares (13 points). The figure shows three 
separate fits for the three subsamples. The black solid line is a fit to the 
gas/stellar dynamics data (black X points), the dashed blue line is a fit to 
the maser modeling data (blue squares) and the dot-dashed red fine is a fit to 
the reverberation mapping data (red triangles). The gas/stellar dynamics fit 
and the reverberation mapping fit ai'e fairly compatible and both close to the 
overall fit. The maser modeling data, however, follows a noticeably shallower 
slope. Given their small sizes it is difficult to be certain if there is really a 
conflict between the subsamples. 

galaxies central bulge, but this is also true for rival theories. 
Thus Pitch Angle, Sersic index, a, and bulge luminosity all 
tend to correlate with each other because they all indirectly 
measure the central bulge mass. There is strong evidence that 
this in turn correlates to central black hole mass. It may also 
ultimately depend on the dark matter halo concentration in 
some way still to be properly elucidated. 

3.3. Subsamples of the Full Direct Sample 

It is instructive to examine the trends in the overall fits of 
the three different measurement techniques we have utilized 
in our data sample. Below, we consider these three techniques 
separately to verify that each individual subsample provides 
consistent results. Additionally, this approach allows us to 
examine any potential differences between active and normal 
galaxies, as two of our subsample groups consist entirely of 
active galaxies. 

Using only the 10 galaxies with mass estimates utilizing 
stellar or gas dynamics, which are mostly in normal galaxies, 
we find a linear fit to the data of 

/og(MBH/M0) = (8.66 ±0.43)- (0.088 ± 0.03 1)P. (5) 

This provides a x^ = 1 -2 and a scatter of 0.39 dex about the lin- 
ear fit (FigureH). This fit differs from the value of lSeigar et alj 
(120081) when only using mass measurements from stellar or 
gas dynamics, but produces a consistent result, though one 
with much higher error, and is similarly consistent with our 
earlier results. 

Examining the masers by themselves (Figure |2]l results in a 
linear fit of 

log{M^n/MQ) = (7.43 ± 0.23) - (0.027 ± 0.012)P. (6) 

The scatter about the fit fine for the masers is 0.30 dex. This 
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data is comprised of 13 galaxies with detectable masers, 12 of 
which have not been measured using stellar and gas dynamics 
methods. 

As we can see from Figure|2l the maser sample has a signifi- 
cantly shallower slope when compared to the direct stellar and 
gas dynamics measurements. This, combined with the good 
quality of the fit and low overall scatter, could imply that a 
separate fit is necessary for this population of active galaxies, 
but it is difficult to say anything for certain with the relatively 
small sample size. 

Turning to the reverberation mapping subsample, (also in 
Figure |2]i we find a linear fit of 

logiM^n/Mo) = (8 .48 ± 0. 15) - (0.072 ± 0.009)P (7) 

with a x^ = 0.93 and a scatter of 0.28 dex for the subsam- 
ple of 12 reverberation mapped masses alone. This is consis- 
tent with the results of Seigar et al. (2008). For an AGN only 
sample of maser modeling data and reverberation mapping we 
find 

logiMYin/Mo) = (8.09 ± 0.45) - (0.058 ± 0.025)P (8) 

with a x^ = 3.68 and a scatter of 0.38 dex for this subsample. 

The stellar and gas dynamics sample has a somewhat 
steeper fit than the overall fit, while the maser modeling sam- 
ple has a somewhat shallower fit. Interestingly, the reverber- 
ation mapping sample splits the difference and ends up very 
close to the overall slope. The steeper fit of the direct, nor- 
mal subsample may be simply attributable to the fact that this 
sample lacks smaller black holes and so it is more difficult to 
determine, on the basis of this sample alone, where the slope 
really lies. 

The reverberation ma pping only result sh ows great similar- 
ity to the original fit by lSeigaret al.l (l2008h and is consistent 
with the fit to stellar and gas dynamics data, with much better 
scatter These are significantly different fits from that found 
using maser data above. They are also consistent with our to- 
tal sample and the resu lts b ased on Mbh - o" measurements as 
we will see in Section 13.41 Meanwhile, the combined maser 
and reverberation mapping sample produces a result consis- 
tent with our other fits, belying any notion that the maser re- 
sults are different because these objects are AGN. 

3.4. Comparing M-P results with M-a 

Here, we consider one of the most common techniques for 
estimating galaxy black hole masses, the M-cr relation. This 
technique is used to estimate the mass of a central SMBH 
by measuring the velocity dispersion of stars in the galactic 
bulge. We consider a set of galaxies using the M-cr relation 
from Ferraresel (12002!) . This data set utilizes a single fit to the 
M-cr relation and allows us to fill in portions of the right hand 
side of the MBH-Pitch Angle relation. We include 23 galaxies 
from this set, 20 of which do not have direct measurements. 
This 20 also includes 3 of the galaxies for which we have mass 
limits based on stellar or gas dynamics, allowing us to replace 
those limits (for the purposes of this section) with the mass 
estimate as derived from the galaxy's a. 

Figure |3] illustrates the results of comparing central veloc- 
ity dispersion (ac) data with spiral arm pitch angle. The left 
panel of Figure[3]demonstrates a tight correlation between cjc 
and pitch angle. The right panel similarly illustrates a corre- 
lation between the 23 masses derived from the M-cr relation 
of iFerrarese, (,2002,) with pitch angle. The fit to the mass vs. 



pitch angle data is 

tog(MBH/M0) = (8.47±O.24)-(O.O89±O.O13)P (9) 

with a x^ = 4.86 and a scatter of 0.48 dex. This result is con- 
sistent with the earlier results based on direct measurements 
of stellar and gas dynamics as well as reverberation mapping 
data. Note that while the scatter is high it is consistent with 
the scatter found in lGiiltekin et al.l (12009) for the total sample, 
^ 0.44 dex, and lower than the result for late-type galaxies 
alone, 0.56 dex. 

Taking the cr-derived masses of the 20 galaxies without di- 
rect measurements and adding them to our direct sample of 
34 gives us 54 galaxies with mass measurements, either direct 
or indirect (Figure 2)- A fit to these data points gives: 

%(Mbh/M0) = (8.36±O.15)-(O.O76±O.OO8)P (10) 



with a x^ = 10.43 and a scatter of 0.4 5 dex (Figure 
result is strikingly close to the result of lSeigar et al.l 
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%(Mbh/M0) = (8.44±O.1O)-(O.O76±O.OO5)P (11) 

To reiterate, the modest change in the M-P correlation re- 
ported in this paper is entirely due to our decision to drop indi- 
rect mass measurements from our sample, because the num- 
ber of available direct measurements has doubled from the 
previous work. Thus our results are consistent with the earlier 
one reported. This, coupled with the fact that the scatter in 
the relation has improved, suggests that the correlation is not 
simply the result of an initially small data set. 

Therefore, we have two fits to the data to which we attach 
special importance. The first is our preferred fit that includes 
all direct measurements and which we adopt as the SMBH 
M-P relation. The second includes all direct data plus val- 
ues of the black hole mass derived from cjc- It is heartening 
that the two fits are consistent with each other, though the 
combined sample has much tighter constraints on the relation 
even with its modestly larger scatter of 0.45, as opposed to 
0.38, dex. 

The fit using only direct measurement data has a shallower 
slope, 0.062 ± 0.009, than the all-data fit, but the two are 
fully consistent. Although we place greater faith in the direct 
data, we must also acknowledge that it is missing a signifi- 
cant amount of data on the right hand side of the relation. Be 
that as it may, we prefer to rely on the fit based only on direct 
measurement data for the final result of our correlation. 

We now examine how consistent the masses generated from 
the M-P relation are with the masses taken from the literature 
(Figure |5]l. Those galaxies with upper and lower limits are 
included in the right panel of the Figure. Notice that M33 is a 
distinct outlier from the rest of the data. 

3.5. Comparisons with Mass Limits 

Finally we look at those galaxies for which only limits are 
available on their black hole masses. We take a look at the 12 
galaxies with upper mass limits set by stellar and gas dynam- 
ics effects, as well as four galaxies which have lower limits 
placed on their masses due to estimates of their luminosity in 
ratio to the Eddington Limit, to check for any inconsistencies 
of the M-f relation with these mass limits. 

As can be seen in the right panel of Figure |5] the limits, 
with few exceptions, are consistent with the resulting fit. In 
the cases where the limits are not consistent, we find that they 
are either still consistent with the scatter observed in the fitted 
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Fig. 3. — Comparisons of pitch angle data with measurements of core stellar velocity dispersion (cr) for galaxies from lFerraresel (120021) . Left: Gc compared with 
our measured pitch angles. Right: Masses taken from the M-ct relation compared to spiral arm pitch angle. Both figures illustrate a strong conelation between 
(a) and spiral arm pitch angle, as one would expect from our argument that both measure the mass of the galaxy's central bulge (see section|A). The fit to the 
SMBH Mass-Pitch angle relation is log{MBH/MQ) = (8.47 ± 0.24) - (0.089 ± 0.013)P with a x" = 4.86 and a scatter of 0.48 dex. 
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Fig. 4. — SMBH Mass-Pitch Angle relation for all available directly mea- 
sured blacli hole masses (as in Figure[T) and for those masse s esti mated indi- 
rectly via cr in our prefen'ed sample(see Table [Tjand Section IzTl for details). 
The fit to the SMBH Mass-Pitch Angle relation for this extended dataset is 
/og(MBH/M0) = (8.36 ±0.15)- (0.076 ±0.008)P with a x^ = 10.43 and a 
scatter of 0.45 dex. This fit is consistent with that obtained from our sample 
of directly measured black hole masses (Figure [T] and is almost identical to 
the fit given in I Seigar et al. 2008). The addition of cr-derived masses, sev- 
eral of which are at the low mass end of our distribution, tend to confirm the 
validity of the complete sample fit shown in Figure [T]against the much shal- 
lower fit found for the maser-modeling-only data shown in Figure|2] Black X 
points represent data from stellar or gas dynamics (10 points), blue squares 
represent data from maser modeling (12 points), red triangles come from re- 
verberation mapping data (12 points) and magenta octagons represent masses 
derived from the M-cr relation (20 points). 

points, or that the pitch angle is at least more consistent with 
the measured a of the galaxy. 

M33 has direct estimates which place an upper limit on its 
mass of 3.0 x 10^ Mq (.Merrittet al.ll2001h . But note that its 
<T suggests a much greater mass than this, one that is more in 
line with our correlation (for illustrative purposes this a de- 
rived value of M33's black hole mass is included in Figured 
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Fig. 5. — Measured black hole mass of our sample of galaxies compared 
to the mass derived by applying the fit illustrated in Figure[r|to our measured 
pitch angles for the same objects. The sohd black fine is included to illustrate 
the 1-1 relation, and the distance of each point from the line gives the residual. 
It will be noticed that some of the limits contradict our fit, as discussed in the 
text. Two in particulai; M33 and IC 342, do so strikingly. 



above). Certainly it would be reasonable to expect, if there 
is a discrepancy between a and more direct black hole mass 
measurements, that the pitch angle would follow a. Similarly 
IC 342 has an upper limit mass, based on gas dynamics, of 
Mbh = 5.0 X 10^ Mq which contradicts what one would ex- 
pect from its measured a (again we included the a derived 
value for IC 342 in Figured above). 

The lower limits placed on the sample using the Eddington 
Limit are also illustrated in the right panel of Figure |5] This 
plot does illustrate that the lower limits placed on the masses 
by their Eddington luminosities are low, but consistent with 
the SMBH M-P Angle relation. 
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3.6. Accuracy ofM-P relation for estimating black hole 
masses 



In iDavis et aP (120 12h we show that pitch angle in many 
cases can be measured to within an accuracy of 3 degrees or 
lower This is especially true for grand design spirals, some of 
which can have pitch angle errors of little over 1 degree. If we 
accept 3 degrees as typical of relatively high quality data then 
the resulting error in black hole mass estimation obviously de- 
pends on the slope of the correlation. A steep slope will trans- 
late a modest pitch angle error into quite a large black hole 
mass error Fortunately the correlation is not especially steep. 
The steepest possible slope consistent with the error in our 
preferred correlation is 0.062 + 0.009 = 0.071. Thus an error 
in pitch angle of 3 degrees translates to a maximum error in 
the log of the black hole mass of 0.071 x 3 = O.213Mbh/M0. 
Thus, if we had a pitch angle which translated into a black 
hole mass of a 1 x 10^ Mq, such an error would mean that the 
mass of the black hole, as estimated by our correlation, could 
vary from 1.2 x lO*" Mq on the high end to 7.87 x 10^ Mq 
on the low end or roughly 1 ± 0.2 x 10^ Mq for the case dealt 
with here; the relative error would be lower for a larger black 
hole. 

In the case of poor quality data or some especially floccu- 
lant galaxies, errors can be 10 degrees or higher In this case 
the error in black hole mass is 0.7 IMbh and so the same black 
hole would have a mass, with error of (1 ±0.71) x lO^'M©. 
Therefore pitch angle measurements with errors greater than 
10 degrees will do little better than estimating black hole 
masses to within an order of magnitude. 

3.7. Sersic Index 



It was shown in iGraham et al.l (1200 Ih that Sersic in- 
dex, which measures the light concentration of a galac- 
tic bulge, correla t es stro ngly with black hole mass and in 
iGraham & Driveri (l2007h a quadratic relation between the log 
of the Sersic index and black hole mass was established. Thus 
the Sersic index provides another observational estimate of 
SMBH masses through images of galactic bulges, and a good 
competitor to the use of Pitch Angle measurements in estimat- 
ing SMBH mass. Certainly, Sersic index is capable of estimat- 
ing SMBH masses in early-type galaxies where Pitch Angle is 
unusable. The Sersic indices for several galaxies in our sam- 
ple, al ong with correspon ding mass estimates from the rela- 
tion of IGraham & Driven (2007) are included in Table |2] We 
will examine the quality of Sersic index-based mass estimates 
using these values as well as the morphologies of the galaxies. 

The left panel of Figure |6] illustrates the relationship be- 
tween Sersic index and the Pitch Angles for galaxies in our 
sample. In the figure the black triangles are spiral galaxies 
without bars, red hexagons are spiral galaxies with weak bars, 
and blue squares are barred spirals. As can be seen from 
the figure there is little apparent correlation between these 
two properties. This is in contrast to the fairly tight rela- 
tion between Sersic index-produced masses compared to the 
masses from the literature, see the right panel of Figure Yet, 
both relations seem to provide reasonable estimates of SMBH 
mass. 

The right panel of Figure |6] compares the masses from lit- 
erature with the masses derived from the Sersic index. The 
same point types are used in the right panel as the left. A 
point of interest seems to be the low scatter about the 1 : 1 line 
for the regular galaxies when compared to barred galaxies, 
but that is not surprising since it is more difficult to measure 



n in barred galaxies using our isophotal fitting technique. The 
scatter for the total sample is ^ 0.60 dex, but is strikingly only 
0.15 dex for spiral galaxies without bars. This would seem to 
suggest that Sersic index is an extremely useful tool in esti- 
mating Mbh for elliptical galaxies and even non-barred early- 
type spirals. For later-type spirals one has to be extremely 
careful and adopt the approach of using a multi-component 
two-dimensional morphological fitting routine (such as GAL- 
FIT dPeng e t al. 2002)). This can become computationally ex- 
pensive, especially when compared to our approach for deter- 
mining spiral pitch angles. 

The results of this work suggest that for a grand design spi- 
ral galaxy the M-P relation, with its low scatter compared to 
other methods when applied to spiral galaxies, could be the 
preferred method of estimating the central black hole mass. 
This may even prove true for all classes of spiral galaxies. 
Even if spectroscopic methods are preferred in some cases, a 
considerable advantage remains that the M-P relation requires 
only imaging data, which is becoming plentiful at high qual- 
ity. 

Accurate Sersic indices still require more image process- 
ing and initial estimates to determine, even when a code such 
as GA LFIT (Peng et al. 2002)) or BUDDA (de Souza et aQ 
I2OO4I) is used for barred spiral galaxies. In these cases, there 
are a larger number of parameters that need to be fit, three pa- 
rameters for the bulge, three for the bar, and two for the disk. 
A large amount of degeneracy may result, making it more 
complex to get an accurate handle on the Sersic index than 
to measure the pitch angle. These methods have bee n used 
to me asure large numbers of galaxies, (Gadotti & KauffmannI 
120091 e.g.), but they still require significant time for image 
processing. As discussed above unresolved nuclear flux will 
create difficulties in determining Sersic index. Automated 
codes may have difficulties disentangling the unresolved nu- 
cleus and the bulge for distant galaxies, this is an issue for 
severa l galaxies in our sample. The method of iDavis et al.l 
( I2OI2I) requires little in the way of image processing in or- 
der to measure the pitch angle of the spiral arms, and requires 
littl e user input . In th e case of late-type non-barred galax- 
ies, JD^ireFal] (HqTII) compare SMBH masses (derived from 
spiral arm pitch angle s using the relation here) with Sersic 
indices measured from IGraham & Driveri (120071). and an en - 
couraging result is found. As discussed in IDavis et al.l (120 13b . 
there are strong reasons to prefer the pitch angle derived mass 
function to the Sersic index derived mass function in the case 
of late-type galaxies. Nevertheless, it is very encouraging that 
the evidence presented here suggests that Sersic index and 
pitch angle estimates of black hole mass are compatible for 
non-barred spirals. This raises the immediate prospect that a 
combination of these two approaches (i.e. using Seric index 
for ellipticals and pitch angles for spirals) could produce a 
black hole mass function for all types of galaxies using imag- 
ing data alone. Additionally, a potential method for rapid 
automate d pitch angle rneasure ments in development is dis- 
cussed in lDavis & HavesI (120121) . 

This comparison provides another opportunity to examine 
our results. By comparing the residuals, derived by subtract- 
ing these results from the masses in literature, we may eluci- 
date whether the correlations are independent. This will also 
inform us as to whether the scatter can be reduced through a 
combination of several parameters. For resi duals derived from 
the lin ear Sersic index mass relation of IGraham & Driveri 
( l2007h versus pitch angle we find a correlation coefficient 
of -0.38 with 80.4% significance. This is a moderate anti- 
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Fig. 6. — Left: Sersic index compared to pitch angle for a subsample of galaxies used in this work. The black triangles represent regular (non-barred) spiral 
galaxies. Red hexagons are galaxies with morphological classifications as weak bars. Finally, blue squares are barred galaxies. Right: A comparison of log(Mass) 
from the literature with the log(Sersic index-Derived Mass) using the relationship of Graham & Driver (2007). The solid line through the data represents the 1:1 
line. Note that non-barred galaxies have much less scatter about the 1:1 fine then either vaiiety of barred galaxy. The scatter for the total sample is ~ 0.60 dex, 
but is only 0.15 dex for spiral galaxies without bars, 0.73 and 0.72 dex for barred and weakly barred galaxies. No Umits are used in this right panel, hence there 
are fewer points than in the panel on the left. 

correlation, but essentially insignificant. Using the preferred 
quadratic relation from Graham & Driver (2007 ), the corre- 
lation coefficient becomes 0.33 (weak correlation) which is 
97.3% significant. This is a greater than 2- cr result. For pitch 
angle residuals versus Sersic index we find a Pearson Rank 
correlation coefficient of 0.74, with a significance of 99.99%, 
about a 4- (T result. 



4. DISCUSSION 

4.1. Comparisons to Other Methods 

In considering the fits to the various categories of data (stel- 
lar and gas dynamics, masers, reverberation mapping, a) it 
is notable that the scatter is less than 0.48 dex in all of the 
samples and subsamples considered. In the case of the di- 
rect measurements using stellar or gas dynamics, the scatter 
is 0.39 dex. This is comparable in scatter to the best of the 
other galactic features which are known to correlate to black 
hole mass. In particular, it is comparable to the scatter for the 
M-cr relation, 0.44 dex, for the full sample of all galaxies in 
iGultekin et al. (2009). The results of Hiiring & Rix (2004) 
produce ^ 0.3 dex scatter using a sample of 30 galaxies. 
These results rely primarily on early type galaxies (Ellipti- 
cals -I- SO) to generate these results, as only 5 of these galaxies 
are spiral galaxies. This provides us with another useful tool, 
as bulge luminosity is another imaging based observable re- 
lation that may be complimentary in a full census of SMBH 
mass. It is true that our correlation holds true only for one cat- 
egory of galaxy (spiral galaxies), but many of the other corre- 
lations hold greater uncertainties and scatter for just this type, 
which is a significant portion of the total population of galax- 
ies. It is further worth noting that the scatter for M-ct, using 
only late-types, is about 0.56 dex, somewhat worse than our 
scatter (Giiltekin et al. 2009). Bulge luminosity and Sersic in- 
dex (both of which depend on measurements of the central 
core of the galaxy) encounter difficulties with spirals, espe- 
cially barred spirals, where one has to subtract the disk and 
bar components to find the true bulge component. Routines 



such as GALFIT and BUDDA provide powerful tools to make 
the necessary measurements of the luminosity from bulges 
JPeng et alj2002l:lde Souza et alj2004l) . These results may be 
use d with relations such as t hose from iHaring & Rix. (.20041) 
and lGraham & Drive J (l2007h to estimate SMBH masses. As 
one can imagine, the more free parameters you have to fit to 
the observed surface brightness profiles (and there are at least 
5 in the case of non-barred spiral galaxies, and as many as 8 
free parameters in barred spirals), the greater the degeneracy 
between these parameters. As a result measuring the value 
of the Sersic index determined from such fits becomes in- 
creasingly complicated. In complex systems (i.e., those with 
bars and disks, compared to pure bulge or elliptical galaxies) 
SMBH mass determinations based on pitch angle also appear 
to be consistent in galaxies with pseudobulges, where other 
relations appear to break down. It is certainly true that many 
of the black hole mass func tions published to d ate concen- 
trate first on early-types (e.g. lMarconi et al1l2004l) . We intend 
to use t he MfiH-pitch angl e relation to redress this imbalance 
(in prep lDavis et al.ll2013l) . Since spirals will typically have 
undergone few major mergers in their history in comparison 
with ellipticals, this technique, which works well for spirals, 
can help develop information which applies particularly to the 
accretion history of black holes, as opposed to the merger his- 
tory. 

The fact that various macroscopic and morphological fea- 
tures of galaxies correlate to each other has been noted since 
galaxies were first observed in any number by Edwin Hub- 
ble in the 1920s. Theoretically, the possibility that the black 
hole at each galaxy's center should correlate with these fea- 
tures is by no means certain but not implausible. A common 
assertion based upon strong observational evidence is that the 
mass of the central black hole correlates to the mass of the 
galaxy's central bulge or core area (e .g Maaorrian et al. 199^ 



Marconi & Huntll2003l: IHaring & Rix.,2004c .Kormendy et aH 
20111) . The reason for such a correlation is not settled, though 



various proposals have been made. One such proposed mech- 
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anism assumes that the depth of the potential in the central 
region in which the black hole resides governs the amount of 
available fuel for accretion. If the black hole becomes too 
large and accretes too rapidly, radiation pressure will force 
fuel out of the well, thus starving itself of further growth. 
Thus, the mass of the central region pl aces an upper limit on 
growth of the black hole (Sil k &Reesl [19981) . As long as the 
SMBH reaches this limit at some point during its life, its mass 
should correlate to the mass of its bulge. Recently, it has been 
suggested that the causal relation is reversed. Instead of the 
mass of the core controlling the growth of the black hole, it 
is the mass of the black h ole which controls the growth of the 
core ( iJahnke et al.l 120091) . Either way, the data suggests that 
such a correlation exists in some form. 

Another observable which seems to correlate well with 
SMBH mass is the Sersic index o f a galaxy (n), as discusse d 
above and in lGraham eTaP (12001): Graham & Driver ('2007]). 
The Sersic index of a galaxy measures the light concentra- 
tion of a galactic bulge, and thus should correlate to th e mass 
concentration of a galaxy. As we saw in Section l3?7l this re- 
lation has potential difficulties when dealing with a large por- 
tion of the population of spiral galaxies. This may be a result 
of the techniques adopted to measure Sersic index in compli- 
cated systems, but this alone may make it difficult to auto- 
mate the measurement of Sersic index, whereas techniques to 
automate the m easurement of spiral p itch angle are already 
being explored ( iDavis & Havesll2bl2h . The relation between 
Sersic index and SMBH mass is however a complimentary re- 
lation that could be used to estimate masses for elliptical and 
SO galaxies as well as confirm results on non-barred bulge- 
dominated spirals. 

Since direct measures of black hole masses are intrinsically 
difficult measurements to make, there is great interest in the 
use of quantities like bulge luminosity, a, and n as markers 
for the study of a supermassive black hole mass function and 
its evolution. There is a level of discomfort with the use of 
such markers that centers around the vast difference between 
the scale of the black hole (which, in terms of the region of 
strongly curved spacetime, is similar in size to our solar sys- 
tem, and in terms of the region within which the black hole 
mass dominates the local stellar mass is only a kiloparsec or 
so) and the scale of the galaxy's central core (on the order 
of 10 kpc). Is it really plausible that reliable correlations be- 
tween quantities on these greatly disparate scales exist? How 
much more cautious should we be in comparing the same 
black hole with a quantity whose scale spans the entirety of 
a disk galaxy, a scale of tens of kiloparsecs? 

Nevertheless, there are excellent observational and theoret- 
ical grounds for believing that the pitch angle of spiral arms 
should correlate well with the mass of the gala xy's central 
bulge despite the significant difference in scale (Lin & Shu' 
1964; Bertinetal. 1989a b; Fuchs 1991; Block etal. 1999; 
Fuchsll2000t ISeigar et al.l 12004 I2005L 120061 12008 ). Indeed, 
not only should this relation be causal, but we can expect it, 
on theoretical grounds, to be quite tight. 

5. CONCLUSIONS 

We have measured spiral arm pitch angles for a sample of 
67 disk galaxies with previously determined SMBH masses. 
The SMBH masses for these galaxies were taken from a va- 
riety of sources which determined the mass via several meth- 
ods including: (1) Direct measurements using stellar or gas 



dynamics; (2) Maser modeling; (3) Reverberation mapping; 

(4) Stellar velocity dispersion and the Mbh - (^ relation; and 

(5) Eddington limits for lower mass limits on the SMBHs. 
The results of several fits to these samples may be found in 
Table|3] Our main conclusions are as follows: 

• Using only galaxies with direct SMBH mass measure- 
ments based upon stellar and gas dynamics in normal 
galaxies as well as maser modeling and reverberation 
mapping in active galaxies, we find a SMBH Mass- 
Pitch angle relation of 

/og(MBH/M0) = (8.21±O.16)-(O.O62±O.OO9)P 

• If we include also select indirect black hole mass esti- 
mates whichwereusedm the previous correlation stud- 
ied in lSeigar etal] (120081) we find 

/og(MBH/Mo) = (8.36 ±0.15)- (0.076 ±0.008)P. 

This rela tion is statis t ically consistent with that pre- 
sented in lSeigaret al.l (120081) . and virtually identical in 
slope. 

• Our scatter is comparable to, or better than, that in the 
M-cr relation. 

• Our technique does not require observationally expen- 
sive spectra. The method is also cosmologically inde- 
pendent, since logarithmic spirals are self similar. 



• Using the relationship of l&aham & Driver! (12007) with 
a sample of our galaxies illustrates that there is more 
scatter in Graham's relationship for barred galaxies. 
Meanwhile, our relationship should not be affected by 
bars. Thus two SMBH mass estimators which can take 
advantage of imaging data only (pitch angle and Sersic 
index) are complementary. 

Equation [TO] is ext r emely useful for measuring SMBH 
masses. In Davi s et alj (120 12h we showed that pitch angle can 
be measured extremely reliably. In future papers we intend 
to use equation [10] to d etermine a l ocal SMBH mass func- 
tion for spiral galaxies (IDavis et al.ll26l3 ), and for higher-z 
spirals, providing we can determine how (or if) the SMBH 
Mass-Pitch Angle relation evolves as a function of look-back 
time. We can often measure pitch angle to within a relative 
error of 3° or less. This translates to a relative error in the 
logarithm of the black hole mass of 4%. 
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Fit Parameters 
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0.38 


34 


1.2 


0.39 


10 


0.9 


0.30 


14 


0.9 


0.28 


14 


3.7 


0.38 


28 


4.9 


0.48 


23 


10.4 


0.45 


54 



Sample Fit X Scatter (dex) Sample Size 

Stellar/Gas, Masers, and Reverberation Mapping'^ /ogCAfBH/M©) = (8.21 ±0.16) -(0.062 ±0.009)P 

Stellar and Gas Dynamics ^^(Mbh/Mq) = (8.66 ±0.43)- (0.088 ± 0.03 1)P 

Masers ^^(Mbh/Mq) = (7.43±0.23)-(0.027±0.012)P 

Reverberation Mapping ^^(Mbh/Mq) = (8.48±0.15)-(0.072±0.009)P 

Masers and Reverberation Mapping logiMBH /Mq ) = (8.09 ± 0.45) - (0.058 ± 0.025)P 

M-CT /og(MBH/MQ) = (8.47±0.24)-(0.089±0.013)P 

AH /og(MBH/MQ) = (8.36±0.15)-(0.076±0.008)P 

Note. — A: Thi,s fit represents our preferred fit to the data. Column 1 Type of galaxy mass measurements used in sample. Column 2 Sample fit. Column 3 x lor 
the fit. Column 4 Scatter about the fit in dex. Column 5 Number of galaxies in the sample. 

APPENDIX 
SPIRAL DENSITY WAVES AND PITCH ANGLE 

In the modal density wave theory it is fairly straigh tforward to see that the pitch angle of the spiral arms must vary inversely 
with the mass of the central bulge of the galaxy. From lBertin & Liiil ( 119961) we have 

m 
tan/=— (Al) 

rk 

where / is the pitch angle of the spiral pattern, m is the number of spiral arms, r is the radial position in the disk, and k is the "local 
radial wave number" which is not constant over the disk but is related to the local wavelength of the density waves by A = In/k. 
Shu ( 1984) deals with the case where the pitch angle is relatively small, which is true for Saturn but also for the tightest spiral 
arms in galaxies, and in the treatment found there one sees how 

1^1 = :r^ (A2) 

where a;, is the surface mass density in the disk and D is an expression which can be understood as the distance, in frequency 
terms, from Lindblad resonance of the gravitational potential. 
A Lindblad resonance occurs when 

uj-mn = ±K (A3) 

where fi is the 9 (or tangential) frequency of the gravitational potential of the central mass of the galaxy (for a point mass this 
is ri = y^GM/P, where M is the mass of the planet and r the radial distance from it) and k is the epicyclic frequency or radial 
frequency of the same potential. For a point mass this is the same as fl. Finally w is a frequency associated with the particle 
orbiting in this potential, in practice likely to be some kind of forcing frequency which encourages the particle to move with a 
frequency not quite that of the main potential. For instance, in the case of spiral density waves in Saturn's rings this co would be 
associated with the perturbing influence of a nearby moon. Since both il and k will vary with radial coordinate r, it follows that 
there will only be certain values of r for which the above relation is satisfied. Such locations are called Lindblad resonances and 
denoted rL. In the theory, spiral density waves emanate from such Lindblad resonances. 
The definition of D is 

D = K^-(uj-mnf (A4) 

and clearly at a Lindblad resonance it follows that D = Q. As a wave moves away from this resonance D will no longer be zero. It 
is easy to show, by a Taylor series expansion, that the key to the pattern of density waves propagating from the resonance is the 
first derivative of D. 

The details of calculating D and its derivative are complex because we are dealing with a system which does not encourage the 
use of the usual simplifying assumptions. The two main frequencies il and k are generated by a decidedly non-point source, the 
central region of the galaxy, and there are many perturbing effects from all of the stars and other material in the disk. Nevertheless 
it is obvious that D (xM where M is the of the central gravitational source (it would be the mass of Saturn in the case of Saturn's 
rings). 

Recall that if we have a point mass source we would have Q = y/GM/r^. In practice, this is far from the case, but we do expect 

that H, = \/GMf{r) where the particular functional dependence on radial distance is uncertain, but where M is known to be all of 
the mass inside the orbit of the particular star or object whose motion we are following. Especially for stars in the inner part of 
the disk, more specifically for density waves emitted from the inner Lindblad resonance, this M will be close to the mass of the 
central bulge of the disk galaxy. 

We must also consider the radial frequency of the orbit (which is, in general, eccentric) and we recall that for an orbit in the 
equatorial plane of an axisymmetric potential 
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In the point mass case K=il, but if il has the generic form given above then we find that 

2 iVGMfir) 



'GM{2rf{r) + r^df{r)/dr) 



■ 2GM 



i^f'^'f'i 



(A6) 



Thus provided il is proportional to VM, we find that k will be also. The frequency uj is not so constrained in general, but near 
the Lindblad resonance it must be close to il and k and so it will also be effectively proportional to VM- Since our density waves 
originate at the Lindblad resonances it follows that this is true in our case. 
We return to the key quantity D and begin from 



D = K^ -(uj-mny = 2GM{2f + rf 



dr - 



Since near resonance we expect cj to have a value such that the second term on the right is small, we rewrite it as 

D = 2GM{2f + rf^) - GMf 
dr 



(A7) 



(A8) 



So we find that D (xM, and the derivative of D with respect to r will also be proportional to M, which is (approximately) constant, 
especially to small changes in radial distance. The fact that D ocM means that 



tan/ oc (Jg/M. 



(A9) 
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